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^.^ Abstract 

o, 

P\J I We analyse dissipation in quantum computation and its destruc- 

^p ■ tive impact on efficiency of quantum algoritlims. Using a general 

Qv . model of decoherence we study the time evolution of a quantum 

register of arbitrary lenght coupled with an environment of arbi- 
o ' trary coherence lenght. We discuss relations between decoherence 

and computational complexity and show that the quantum factori- 
C I sation algorithm must be modified in order to be regarded as efficient 

^ ' and realistic. 



Jx^ '• 1 Introduction 



Quantum computers can accept input states which represent a coherent su- 
perposition of many different possible inputs and subsequently evolve them 
into a corresponding superposition of outputs. Computation, i.e. a sequence 
of unitary transformations, affects simultaneously each element of the su- 
perposition generating a massive parallel data processing albeit within one 
piece of quantum hardware. As the result quantum computers can efficiently 
solve some problems which are believed to be intractable on any classical 
computer (Deutsch 1985, Deutsch and Jozsa 1992, Bernstein and Vazirani 
1993, Simon 1994, Shor 1994). The most striking example is the factoring 
problem: to factor a number N oi L digits on any classical computer re- 
quires an execution time that grows exponentially with L (approximately 
expL^/^ for the best known algorithms such as the Number Field Sieve 
(Lenstra et al. 1990)). In contrast Shor (1994) has shown that quantum 
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computers require an execution time that grows only as a polynomial func- 
tion of L (^ L^). A noteworthy consequence for cryptology is the possibility 
of breaking public key cryptosystems such as RSA (Rivest et al. 1979). In 
this context a practical implementation of quantum computation is a most 
important issue. 

This paper shows that Shor's algorithm, in order to be efficient from the 
computational complexity perspective, must be modified (i.e. supplemented 
with an exponentially efficient quantum error correction). This is because 
the error rate, which is due to the interaction with the environment (a ther- 
mal reservoir), grows as an exponential function of L regardless the type 
and the strength of the computer - environment coupling. Consequently, 
in order to obtain factors of N with a prescribed probability of success one 
has to repeat the computation k times, where k grows exponentially with L. 
Let us stress that this observation is not related to any technological limi- 
tations but is a consequence of the computational complexity requirements 
imposed on any realistic model of computation. Any technological progress 
in suppressing dissipation will, of course, help to factor bigger and bigger 
numbers but will not change an exponential increase of the error rate. 

Our analysis is quite general and covers cases where qubits interact with a 
single reservoir and the distance between the qubits is comparable with the 
correlation length of the reservoir. In the limit of a short correlation length 
(i.e. coupling to the independent reservoirs) our results coincide with those 
by Unruh (1995). Also we provide a plausibility argument that any model 
of dissipation will lead to errors that grow exponentially with L. 

The structure of the paper is as follows: Section 2 links the computational 
complexity of randomized algorithms with dissipation. Section 3 outlines ba- 
sic features of the dissipation model we adopted for this paper (the model 
is known under the label "decoherence" ) , Section 4 explains details of de- 
coherence of a single qubit using both semi-classical and fully quantised 
approaches. Section 5 extends the fully quantum analysis to the two qubit 
case first and then to the general case of an L-qubit register. The results 
are summarised and discussed in the last section. 



2 Errors in Randomized Algorithms 

An algorithm is said to be efficient if its execution time increases no faster 
than a polynomial function of the size of the input. For example, any 
efficient factoring algorithm will factor a number A^ of size L = log A^ (the 
base of the logarithm is 10 when the size is measured in decimal digits and 
2 when it is measured in binary digits) in time proportional to a polynomial 
function of L. Algorithms do not have to be deterministic; random decisions 
can be taken at certain junctures in the execution of the algorithm. The so- 



called randomized algorithms may sometimes produce an incorrect solution, 
however; they have the property that the probability of error can be made 
arbitrarily small. 

Consider, for example, an efficient factoring algorithm which runs success- 
fully only with probability 1 — e (0 < e < 1). Clearly we know when the 
algorithm is successful because we can check it efficiently by a trial division. 
Assuming that e is independent of the input the success probability of the 
algorithm can be amplified (boosted) arbitrarily close to 1 by repeating the 
computation several times. After k runs the probability of having at least 
one success is 

P{k,e) = l-e\ (1) 

If we fix P{k, e) (i.e. when we agree upon the desired probability of success- 
ful factorisation) then the number of runs k depends on the error rate e. In 
order to retain efficiency of the algorithm we may let e increase with L but 
not faster than a polynomial. Indeed, Shor's quantum factoring algorithm is 
of this type; in a single run it provides factors of A^ with probability 1 — e(L) 
and the probability of error e{L) grows only as a polynomial function of L. 
However, errors in Shor's algorithm are of purely "mathematical" origin and 
correspond to the random decisions which rely on the prime number theo- 
rem (i.e. on the distribution of prime numbers). In a realistic scenario the 
probability of error e{L) has also contributions which are due to the com- 
puter - environment coupling. From the computational complexity point of 
view it is irrelevant how weak this coupling is for a fixed L, it is its depen- 
dence on L which matters. If the interaction with the environment causes 
an exponential increase of the error rate i.e. e{L) = 1 — Aexp{—aL), where 
A and a are positive constants, then k{L) ^ exp(aL) and the randomized 
algorithm cannot be regarded as efficient any more! 

Indeed, our results show that the computer-environment interaction leads 
to this unwelcome exponential increase of the error rate. Shor's algorithm, 
as it stands at the moment, is not robust to the environmental noise and 
requires modifications. This statement, however, should not be taken as a 
fundamental argument against possibilities of experimental implementation 
of quantum factorisation. It simply shows that further efforts, both theo- 
retical (e.g. designing an exponentially efficient quantum error correction) 
and experimental (e.g. customizing the environmental noise), are necessary 
to make quantum computation practical. 



3 Decoherence 

Let us consider a quantum register composed of L qubits; each qubit being 
a two-state quantum system described by a two-dimensional Hilbert space 



7^2 with computational basis states labeled as | 0) and 1 1). The 2^ dimen- 
sional Hilbert space of the whole register is an L term tensor product of 
0,2 pertaining to each qubit. Any quantum state of the register can be 
described by a density operator of the form 

pit)= Ep^.WN)0'I' (2) 

where the computational basis | i) is defined as a tensor product of the qubit 
basis states: 

I i) = I iL-i) ® I iL-2) ® . . . ® I ?o) • (3) 

The r.h.s. is the binary decomposition of number i = Y.t^o 2'i;, e.g. for a 
two qubit register | 2) = 1 1) | 0). 

Quantum computation takes its power from quantum interference and en- 
tanglement. The degree of the interference and entanglement in an L qubit 
register is quantified by the coherences i.e. the off-diagonal elements of the 
density operator in the computational basis pij {i 7^ j). When a quantum 
computer is in contact with an environment, here modeled as a thermal 
reservoir, the resulting dissipation destroys the coherences and changes the 
populations (the diagonal elements). When a reservoir is in a thermal state 
the coherences change on a much shorter time scale than populations and 
it is decoherence which is responsible for destroying quantum parallelism. 
However, it should be mentioned that in some cases populations can be more 
fragile than coherences. The time scales depend on the type of the system- 
environment coupling and on the quantum state of the reservoir (which 
does not have to be thermal), for example, in some optically driven systems 
coherences change more slowly than populations (Kimble & Mabuchi 1995). 

The notion of decoherence can be intuitively understood as follows. Con- 
sider a computer and a reservoir, initially both in pure states: 

computer 
/- " 1^ reservoir 



l*(0))=(^Ec.(0)N)j® ko) . (4) 

A unitary evolution of the composed system results in an entangled computer- 
reservoir state which can be written as 

I^W) = Ec.WI^)®l«.(t)), (5) 



where, in general, (ttj | Uj) 7^ for i 7^ j. The elements of the register 
density matrix evolve as 

p,,(0) = q(0)c*(0) -^ pij{t) = c,{t)c*{t) {a,{t) I a,{t)). (6) 

In a popular model of decoherence (Zurek 1991), where the environment 
effectively acts as a measuring apparatus, the evolution affects the reservoir 
states (I ftj)} which become more and more orthogonal to each other whilst 
the coefficients {q} remain unchanged. Consequently the off-diagonal ele- 
ments pij disappear due to the {aiit) \ aj(t)) factors. We will analyse this 
type of decoherence in detail in the following sections. 

The characteristic time for the off-diagonal elements to disappear is known 
as the decoherence time. Its value depends on the type of qubits and their 
interaction with the environment and can vary from 10'^ s for nuclear spins 
in a paramagnetic atom to 10~^^ s for electron-hole excitations in a bulk 
of a semiconductor (for a table of decoherence times for various physicsl 
systems see for example DiVincenzo 1995). 



4 Single Qubit Dephasing Mechanism 

In this section we want to describe in detail the basic features of the de- 
coherence model we adopted for this paper. For this purpose it suffices to 
analyse an interaction between a single qubit and the reservoir. The case 
involving two and more qubits will be analysed in the subsequent sections. 

The density operator for a single qubit has the form 

pW= Ep^.WI^)0'I. (7) 

i,j=0 

and we are interested in the time evolution of the coherence poi{t) = p{o{t). 

4.1 Semiclassical Picture 

In a semiclassical analysis of the single qubit decoherence it is assumed that 
the reservoir is a fluctuating classical magnetic field. This approach excludes 
many interesting quantum features, most notably the qubit-environment 
entanglement; however, we have decided to include it here because it can 
simplify some calculations. 

We describe our physical qubit as a fictitious spin 1/2 system following 
the standard correspondence between the density operator p and the spin 



vector s, also referred to as the Bloch vector (Allen & Eberly 1975, Colien- 
Tannoudji et al 1977) 

Pit) = hi + s{t)-^) = la+ E 5«(^K), (8) 

■^ ■^ a=x,y,z 

where (Tq, a = x,y,z are the Pauh spin matrices. Components s^ and Sy 
correspond respectively to the real and imaginary parts of poi) and s^, the 
so called population inversion, is the difference between the two diagonal 
terms pn — poo- We are mainly interested in Sx and Sy because they describe 
the internal coherence of the qubit. 

Consider the time evolution of s in the presence of an external magnetic field 
Bq, which points in the z direction. We label the two eigenstates of the az 
spin operator with |0) (spin down) and |1) (spin up). These eigenstates are 
separated by an energy gap ujq = qBq, where g is the gyromagnetic ratio. 
Hence the spin vector rotates along the z axis with an angular frequency ujq 
(see Fig. 1). 

Let us suppose that the external noise can be modeled as a stochastic time 
dependent field Bz{t). Then the xy component of the spin vector, i.e. Sxex + 
SyCy, rotates around the z axis with an instantaneous angular frequency uJo + 
fi{t), where fi(t) = gB-zit). After time t the awerag'e phase factor acquired by 
Sx and Sy due to the presence of external noise becomes (exp \i j^ Vt{t')dt' \ 
where the average is to be taken over all the possible field configurations. 

In order to evaluate the average over the field configurations it is convenient 
to decompose the stochastic field into a Fourier series as follows: 

Bzit) = Y, ^ke*^'* + ^kC"^^** (9) 



where 6k, &k ^^^ Gaussian time independent random variables. This allows 
us to write 

exp h r Q{t')dA\ = H /exp (-ig H h^e''^'^'' dt' + c.c.|\ (10) 

which shows how the overall phase factor can be expressed as a product 
of independent contributions from each of the Fourier component of the 
fluctuating field. These contributions add incoherently and destroy the spin 
coherences. 

The process of decoherence described above is similar to what is knonw as 
T2 decay in nuclear magnetic resonance (Allen & Eberly 1975). This is to 
be distinguished from Ti decay, which implies a damping of the spin energy 
into the field. Our model does not describe Ti processes. 
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To illustrate the destruction of coherences in the semiclassical picture we 
simulate the energy fluctuations by explicitly introducing the random energy 
fluctuations Bz(t) while calculating the time evolution of s. We do this 
numerically in discrete time steps At. For each time step we select a random 
number k between and 1. If k < 0.1 then E^(t + At) = B^^t) + B', 
where B' is some constant fleld. Correspondingly, if k > 0.9 then Bz{t + 
At) = Bz{t) — B' . Otherwise B remains unchanged. We scale our variables 
dimensionless, set B' = O.lBg and set 5^ = 1. In flgure 2 we show examples 
of such simulations. In a more realistic model the size of fluctuations should 
be randomized as well, but this simple model with flxed-size fluctuations is 
good enough for our illustration purposes. 

Each stochastic calculation of time evolution will produce a practically 
unique history for s. These histories form an ensemble, and if we aver- 
age over its all possible members, we obtain a statistical average which then 
corresponds to the one shown in (|10|). In flgure 2(a) and 2(b) we show two 
examples of single ensemble members, and in flgure 2(c) we show an en- 
semble average of 100 such random ensemble members. We see that as the 
fluctuations average out, the changes introduced to the coherences Sx and 
Sy remain. And as seen in flgure 2(d), with 500 members in the ensemble, 
they eventually disappear altogether as t becomes large, i.e., t ^ Uq^. 

Our semiclassical simulations for a single quantum bit demonstrate how the 
full statistical result is obtained by adding the individual members of the 
ensemble. It should be noted that in the time evolution of a single member 
of the ensemble the amplitude for the oscillations in Sx and Sy remains 
constant. Hence there is no real decay of the coherences present in a single 
ensemble member, but merely a shift in phase. 

In the semiclassical picture it is possible to regard a single member of the 
ensemble as a fully physical entity. It could describe e.g. collisionally in- 
duced energy fluctuations in qubits formed by ions in a trap; these arise 
when the atoms of the background gas collide with the ions and perturb the 
ionic energy level structure. However, when we consider the effect of the 
thermal and vacuum fluctuations, such a semiclassical view on the interac- 
tion of the qubit with the reservoir is not adequate, and one has to apply a 
fully quantum approach, which we shall present next. 

4.2 Fully Quantised Model 

The above description of our decoherence mechanism can be rigorously for- 
mulated in quantum terms. The quantised environment is modeled as a 
continuum of fleld modes. The state of the combined system i.e. qubit -|- 
environment is described by a density operator g{t) which at time t = is 



assumed to be in the form 

^(o)=p(o)®n^kT, (11) 

k 

where R^t is the thermal density matrix of the k mode of the field and 
we have taken advantage of the fact that the density operator of the field 
in thermal equilibrium factorizes into the tensor product of the density 
operator of each of its modes. 

We assume that the dynamics of a single qubit + environment is described 
by the following Hamiltonian 

H = \a,ujQ + Y. h^^^ujk + E ^-(^k&I + ^k^k) (12) 



where the first and the second term on the r.h.s. describe respectively the 
free evolution of the qubit and the environment, and the third term describes 
the interaction between the two (6k; 6J^ are now creation and annihilation 
field operators, here and in the following we put ^ = 1). Hamiltonian (12) is 



equivalent to the one introduced in connection with the tunnelling problem 
(Leggett et al. 1987, Gardiner 1991) and used to model decoherence in 
quantum computers (Unruh 1995). 

Since [a^, H] = the populations of the qubit density matrix, p{t) = Ti jig(t) 
are not affected by the environment and the coupling with the environment 
in our model simply erodes quantum coherence. This means that also in 
our fully quantised model there is no exchange of energy between qubit and 
environment and consequently no Ti type of decay takes place. This is not 
a limitation of the model since, as will became clear in the following, it 
describes adequately the decoherence mechanism which usually takes place 
on a shorter timescale than the energy dissipation. Furthermore this model 
is exactly soluble and allows a clear analysis of the mechanism of entangle- 
ment between qubit and environment which is believed to be at the core of 
most decoherence processes. 

We will postpone a detailed discussion on the form of the coupling g^ be- 
tween qubit and field modes which will depend on the specific characteristics 
of the physical system. 

In order to study the time evolution of our system it will be convenient to 
move to the interaction picture, where the time evolution operator takes 
the form 



U{t) = exp|-z|^*^a,((7k6Le*"'=*' + ^kM-'"'*')rft'| 



(13) 



exp a4E(feW-^kC(t); 



with 

Ut) = 2^k . (14) 

Here U{t) can be described as a conditional displacement operator for the 
field, the sign of the displacement being dependent on the logical value of 
the qubit. In particular for any pure state | \1') of the field 



k 

f/(t) I 1) ® I il/) = |l)®n^(+|aW)|^) (15) 

k 

where the displacement operator -D(^k) is defined as 

Di^^)=exp{bU^-b^a}. (16) 

The above discussion makes it clear that U(t) induces entanglement between 
qubit states and field states. For example, if at time t = the state of the 
system composed by the qubit and the k mode of the field is the tensor 
product of a general qubit state times the vacuum state for the field mode, 
then the interaction Hamiltonian will, at time t, generate an entangled state 

(co|0) + ci|l)) ® |0k) ^ Co|0)| - lUt)) + ci|l)| + lUm (17) 

where \^C,\i(t)) is a coherent state of amplitude |Ck(^)- Such correlations 
between qubit and environment cannot be accounted for in the semiclassi- 
cal approach outlined at the beginning of this section. It is precisely this 
entanglement which, in the fully quantised description we have adopted, is 
responsible for the decoherence process. Indeed the off-diagonal elements 
of the reduced density matrix of the qubit decay due to the fact that the 
overlap between the different field states with which the qubit becomes en- 
tangled diminishes in time. We can formulate in rigorous terms the analysis 
we have outlined above by taking into account all the field modes in ther- 
mal equilibrium. The matrix elements of the reduced density operator of 
the qubit are defined as 

p,,{t) = {z\TTnU{t)g{0)U'\t)\j) . (18) 



Using Eg. (p!5|) we check that, as anticipated, poo(t) = poo(O), Pii(t) = pii(O) 
and for the coherence pio we obtain 



k 

= e-rWpio(O). (19) 

Equation ([T9|) is exact, i.e. no approximation has been made to obtain it. 

To evaluate exp{— r(t)} we have to calculate the average value of the dis- 
placement operator for each mode in a thermal state. This is also known as 
the symmetric order generating function for a harmonic oscillator in thermal 
equilibrium (Hillery et al. 1984, Gardiner 1991). It can be shown that 

rrk{i?kTl^(a)} = exp |-M coth (^) I , (20) 

(we put the Boltzmann constant ks = 1). Thus, in the continuum limit 



r(t) oc / (iklfifb-Pcoth ( — - I ^ 

^ ' J '^ ' \2T) ujI 



oc duj—G{uj)\g{uj)Wl + 2{n{uj))T) ^^ (21) 



.^G(^)|,(^)P(l + 2(n(^)ML^^-^ 



where G{uj) is the density of modes at frequency ui (we have dropped the 
superfluous index fc), {n{uj))T = exp(— ti;/2T)cosech(a;/2T) is the average 
number of field excitations at temperature T and {dk/du) is the dispersion 
relation. In Eq.(^) we can separate the effects due to thermal noise from 
the one due to purely quantum fluctuations. This formal separation allows 
us to identify the existence of various timescales in the decoherence process 
(see also Unruh 1995). Let us first note that thermal fluctuations can affect 
the qubit dynamics only for times longer than the characteristic thermal 
frequency T. For t < T^^ only vacuum quantum fluctuations contribute to 
the dephasing process. 

Furthermore the quantity G{uj)\g{uj)\'^ is in general characterized by a cut- 
off frequency whose specific value depends on the particular nature of the 
physical qubit under investigation. For example if the noise field is a phonon 
field the natural cutoff can be identified with the Debye frequency. More 
generally we can think of the cutoff as due to some characteristic length 
scale in our system below which the qubit-environment coupling decreases 
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rapidly. Therefore we assume {dk/dLu)G{uj)\g{uj)\'^ oc cu^e"'^/'^''. The expo- 
nent n will depend on the number of dimensions of the field. In our analysis 
we will concentrate our attention to the case of one-dimensional field, for 
which n = 1 and of three-dimensional field, for which n = 3. It is evi- 
dent that the so-called quantum vacuum fluctuations will contribute to the 
dephasing process only for times t > uj^^. 

We can identify three time regimes of decoherence: 

• a "quiet" regime, for t < uj~^, where the fluctuations are ineffective in 
the decoherence process 

• a quantum regime, for u^^ < t < T^^, where the main cause of 
coherence loss are the quantum vacuum fluctuations 

• a thermal regime, for t > T^^, where thermal fluctuations play the 
major role in eroding the qubit coherence. 

In order to have some semiquantitative picture of the time dependence of 
the decoherence we need to specify the frequency dependence of the density 
of states and of the coupling. The case of one-dimensional field (n = 1) has 
already received some attention in literature, in this case 

r(t) oc / rf.;e-/- coth (^) l^^^^M. (22) 

An analytic solution of (|22|) can be obtained in one dimension in the low 
temperature limit {uc ^ T): 



T{t) ocln(l + <r) + 21n 



sinh(7rTt) 



■nTt 



(23) 



The first term arises from the quantum vacuum fluctuations while the sec- 
ond is due to the thermal ones. Expression (|23|) reduces to r(t) ~ ult^ for 
t < uj-^, V{t) ~ 2 \nuct for uj-^ < t < T-\ and T{t) ~ Tt for t > T~\ 

These three regimes can be easily identified in Fig.(^), which shows the 
decoherence of a single qubit induced by a one-dimensional field for the 
particular choice uJc/T = 100. 

It is also interesting to consider the three dimensional field case (n = 3), 
where the integral can be evaluated exactly: 



r(t) oc f du ue-'^^'^^ coth (^^ [1 - cos(c<;t)] 
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oc 2r2 2C 



+UJ' 



OJa 



T 

2,— (l + iwct) 

UJr. 



2, — (1 — iojct) 



+ 



(1 + iuJctY (1 - iuJctY 



(24) 



where ({x,y) is the generahzed Riemann zeta function. In Fig.@ is shown 
the decay of a single qubit as a function of time and as a function of the 
ratio 1] = uJc/T in the case of a one-dimensional (a) and three-dimensional 
(b) fields. In (a) the coherence always decays to zero. For r^ ~ 1 this decay 
is dominated by the thermal fluctuations of the reservoir and it is therefore 
exponential. However, for large r] there is an intermediate region where 
the vacuum fluctuations dominate and the decay is roughly linear before 
the thermal regime takes over. In all cases there is a short t region where 
almost no decay is present. As t] decreases, the decay onset moves towards 
larger t. This is because the cut-off frequency ujc determines the extent of 
the "quiet" regime. In (b) the system shows a very different behaviour: 
the decoherence saturates to a value determined by rj. This difference in 
behaviour is reminiscent of the dimensional differences observed in various 
models of phase transitions (Ma 1976) and is due to the suppressed influence 
of low frequency fluctuations in three dimensions. It should however be 
pointed out that at longer timescales, when Ti processes involving exchange 
of energy between qubit and environment take place, this residual coherence 
will disappear. 



5 Decoherence of Quantum Registers 

In the present literature on decoherence processes in quantum computation 
it is usually assumed that in a register of length L each qubit interacts 
individually with a different reservoir. In this case all the analysis of the 
relevance of dissipation from the complexity viewpoint is done assuming 
that Tiit) ~ LV{t) (see next section). In this section we discuss in which 
circumstances this assumption is justifled and we analyse the consequences 
of collective interaction on the complexity analysis, when such effects need 
to be taken into account. 

We will start by considering a system of two qubits at positions r^, r^ re- 
spectively. Their dynamics is described by the following Hamiltonian 



H 



+ 



iaX + K^o + E^L&k^'^ 



k 



<{9lhl + gl*h^) + al{gih 



9'a: 



(25) 
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Hamiltonian (^) is a straightforward generalization of (12) where 5'k,5'k 
are position dependent couphngs. In ( ^5]) we are not considering any direct 
interaction between the two qubits which would be necessary for the condi- 
tional dynamics (Barenco et al 1995). Here we concentrate our attention to 
the simplest case of a system of two qubits coupled only to the environment. 

In the interaction picture the time evolution operator takes the form 

U{t) = exp 1 1 ^ 6t [^a^.(^) + ^6^5 (^)] _ ^^ J^a^a*(^) ^ ^5^5* (^)] 1 (26) 



Here again the unitary evolution produces entanglement between register 
states and environment states. It should be noted however that U{t) acts as 
a displacement operator on the field with a displacement amplitude depend- 
ing on the logic value of both qubits of the register. The following simple 
example will help to clarify this point. Let us consider the following two 
initial states 



|<^^"^) = (cio|la,06)+coi|Oa,U))®|Ok), (27) 

|<^^+^) = (Coo|Oa,06)+Cn|la,U))®|Ok), (28) 



where we consider for simplicity only the vacuum state of the k mode and 
where the Cij are arbitrary complex amplitudes. Due to the qubit-field 
interaction the system will evolve into the states 



i<i>(-)(t)) = cio|ia,o,)i + |(e^-e^)) + coi|o„u)i|(ee-e^)), (29) 

l^^+n^)) = Coo|0„0,)|-i(e^ + ^^)) + Cn|la,U)|i(e^ + e^)). (30) 



Since the amount of decoherence is measured by the overlap between the 
two different field states with which the qubit states become entangled, the 
states I $*•''''), 1$ ) will be characterized by different decoherence times. In 
particular, in the limit r^ k, n, we have ^^ ~ ^^ and therefore in the "sub- 
decoherent" state l"^^"-*) the qubits are not entangled with the field, while 
in the "superdecoherent" state |$*^+)) the qubits are collectively entangled 
with a coherent state with an amplitude which is twice that of the single 
qubit case. 

We can formulate in more rigorous terms the analysis we have just sketched 
above by taking into account all the modes of the field in thermal equi- 
librium. We will again concentrate our attention on the reduced density 
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matrix of the two qubit system whose matrix elements can be conveniently 
expressed as 

Piaja,idM = {ia,ib\Trii{g{t)}\ja,jb)- (31) 

It is straightforward to verify that the analogous of Eq. ([19D is 

k 

Matrix elements of the form Pi^i^^i^j^ show no collective decay, in fact 

k 

which for % = jb reduces to pi,i,,i,i,it) = A,i„,i,i,(0) while for if, ^ jh we have 
Piaia,ibjbi^) ~ Piaia,ibjbi^)^~^^^^ ■ Aualogous results hold of course for matrix 
elements of the form Pi^j^^i^i^- 

Collective decay is instead apparent in the matrix elements of the form 
Piaja,ibjbit) with ia i- ja, ib 7^ jb ■■ 

PloMt) = Pl0,10(0)n^^k{i?kT/^(ek + ek)} 

k 

= e-r^Wpio,io(0), (34) 

Pio,oi(t) = Pio,oi(0)n^^k{i?kT/^(ek-ek)} 

k 

= e-r"Wpio,oi(0). (35) 

The expression for F^ can obtained from ( pOD with the obvious substitution 
{k ^ Ck i ^k- If W6 suppose g^ = 5'ke*^'"" and g^^ = 5'ke*^'"*' , then, in the 
continuum limit, the collective decay rates are 



r±(R,t) oc /rfkl^krcoth (^) L^^fi ± eos(k • R)] (36) 



where R = r^ — r;,. 



The origin of the position dependence of the collective decay constants 
r^(R, t) can be easily explained both in semiclassical and in quantum terms. 
We have already mentioned that in a system of two qubits the entanglement 
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with the field will in general depend on the qubit separation. In a semiclas- 
sical picture if the distance between the two qubits is smaller than the cor- 
relation length of the bath they "feel" the same fiuctuation. More generally 
the random phase-shifts due to the bath fiuctuations of wavelength shorter 
than the inverse of the qubit separation are the same for both qubits. This 
produces a reduction - or an enhancement - of the decay rates depending 
on whether the phase shifts on the two qubits add or subtract. 

In order to avoid possible confusions it must be stressed that the phe- 
nomenon of superdecoherence vs. subdecoherence that we have just de- 
scribed is not the same as the process of superradiance vs. subradiance 
more commonly encountered in litterature (see e.g Allen & Eberly 1975). 
While superdecoherence is due to collective entanglement between qubits 
and environment with no exchange of energy, superradiance is a process of 
collective radiation by a group of closely spaced atoms. Indeed the condi- 
tions for the two processes to occur are quite different: superdecoherence 
requires the qubit to be in a region of space smaller than the coherence 
length of the environment, which is a much more stringent condition that 
the one needed for superradiance to take place, i.e. that the atoms are at a 
distance smaller than the wavelength of the resonant field modes. Further- 
more, it must be noted that some superradiant states are subdecoherent. 
In general the set of states which are robust against the dephasing action 
of the environment will depend on the specific form of qubit-environment 
coupling. 

In the one-dimensional and three-dimensional field cases, the collective de- 
cay constants are, respectively, given by 

^ I ,^ , ^ f°° , ,,/,, , / u \ 1 — cos out , 

sin{ujts 



TfJR, t) oc 2 / rfcue-"/"^ coth ( — ) ^[1 ± cosMJ] (37) 

JO \2T J CO 



r^o(R,t) oc 2 r duiue-'^/'^^' coth (^^ (l-coscut) [ 1 ± 



ut,. 



(3J 



where we have introduced the transit time tg {ootg = k ■ R). 

The analysis carried out for a two-qubit register can be easily extended to 
the case of a register of L qubits, whose density matrix elements can be 
written as 



P{in,jn}it) = {ii-i,iL-2-;io\TrR{g{t)}\JL-i,JL-2,-;jo) (39) 

= PK,i„}(0) 

X 



k 

••••^[(^o-Jo)ef}. 
15 



From the viewpoint of the following complexity analysis it is interesting to 
consider the limiting case in which all the qubits of the register are in the 
same position. In this case all the ^^ are equal and the matrix elements 
with the fastest decay are /5{i„,o„} and P{o„,i„} for which we have 



P{lnfin} =P{i„,o„}(0)' 



,-L^r{t) 



(40) 



In general, 



P{in,j„}W = P{in,jn}iO)Y[Trk\RkTD 



PKjn}(0)exp 



L-1 



X-1 



J2 (^" ~ ^")^k 



.n=0 



n=0 



J2 i^n - in) 



r(t) 



(41) 



which should be compared with the expression for the decoherence in the 
case of independent reservoirs 



L-l 



PKJn} W = PKJn}(0) exp <^ - ^ |(i„ - jn)\T{t) 

I n=0 



(42) 



where '}2n=l \i}n — jn)\ is the Hamming distance between the two qubit 
states. In the independent reservoir case therefore in the worst case the 
density matrix elements decay as exp{— Lr(t)}. 



6 Discussion and Conclusions 



Our analysis shows that decoherence destroys quantum interference and en- 
tanglement in quantum computers, thus decreasing the probability of suc- 
cessful computation exponentially with the input size L. Reservoirs with 
large coherence lengths introduce only asymmetry in the decoherence rates 
of various off-diagonal elements of the density matrix with the same Ham- 
ming distance between its indices (compare Eq. (^Tj) with Eq. (|^)). In 
this case while some coherences are not affected by the decoherence some 
others are destroyed more rapidly. Assuming that all coherences are equally 
important and that for the complexity analysis we choose the worst case 
(i.e. the off-diagonal element which decays at the fastest rate) then we 
are forced to conclude that reservoirs with a large coherence length can be 
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more of a hindrance than a help. In the independent reservoirs hmit the 
fastest decoherence is proportional to exp(— Lr(t)) whereas for the reser- 
voirs with a large coherence length there are off-diagonal elements which 
decay as exp{—poly{L)T (t)) . However, if one can support quantum compu- 
tation with some selected coherences then a large coherence length can be 
utilised. 

Let us suppose that we can manufacture in our laboratory a quantum reg- 
ister of 2L qubits composed of pairs of qubits close enough to each other so 
that each pair is effectively interacting with the same reservoir. Different 
pairs can interact with different reservoirs although the results we are going 
to illustrate are not modified if all the qubits interact with the same reser- 
voir. The idea is that if we use a pair of qubits to encode a bit we might 
effectively decouple the register from its environment. To give a simplified 
example of how this might work let us consider a register of 4 qubits in 
which the two pairs of qubits interact with two different reservoirs. The 
interaction will evolve the system as 

where for simplicity we have considered only the vacuum state of one mode 
for each reservoir, labelled a, b, respectively. It is obvious that states for 
which ia 7^ ja and ib 7^ jb do not get entangled with the reservoir and there- 
fore they mantain their coherence. We might therefore use the following 
encoding 



|6) = |o,i), 

|i) = |l,0). (44) 



This " noiseless" encoding does not pose complexity problems since to encode 
L qubits we need simply 2L qubits. Of course, the experimental realization 
of a register with the above characteristics encounters practical problems; 
it might, for instance, be difficult to excite and manipulate only the stable 
coherences. The main point we want to make, however, is that appropri- 
ate encoding could be a way to overcome some of the problems posed by 
couplings with external noise. 

Is the computational complexity analysis dependent on the character of de- 
coherence i.e. on the type of coupling with the environment? We believe 
that as long as the reservoir is in a thermal state, any model of decoherence 
will lead to the same conclusions. In the independent reservoirs limit the 
fastest decoherence is proportional to exp(— Lr(t)) regardles the functional 
dependence r(t). For reservoirs with a large coherent length the fastest 
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decoherence will be proportional to exp{—poly{L)r'{t)) leading to the same 
unwelcome exponential increase of the error rate. However, let us emphasise 
again that this should not be taken as a fundamental argument against pos- 
sibilities of experimental implementation of quantum factorisation. Quan- 
tum error correction, modification of the existing quantum algorithms, cus- 
tomizing the environmental noise, and probably many other techniques can 
be employed to improve the efficiency of quantum algorithms. 

Even if quantum algorithms for factorisation are not efficient from the com- 
plexity theory point of view, they may still permit factoring of numbers 
which are much bigger than those that can be factored using classical algo- 
rithms. Technological progress in isolating quantum computers from the en- 
vironment and suppressing decoherence will increase the size of the biggest 
number that can be factored by such computers. For if the elementary com- 
putational step takes time r and t is a decoherence time of a single qubit 
then the requirement for the coherent computation to be completed within 
the decoherence time of the computer can be written as 

rL^ < t/L. (45) 

The l.h.s. represents the time needed to factor number of size L with the 
Shor algorithm and the r.h.s. represents the decoherence time of L qubits 
(in the independent reservoirs approximation). From the equation above 
follows that the size L is bounded by 

L<(-) . (46) 

Thus the ratio t/r, which depends on the technology employed, determines 
the limits of the algorithm and it is unrealistic to assume that this ratio can 
be made infinite. For a review of possible values of t and r see DiVincenzo 
1995. 

Clearly the decoherence problem calls for a satisfactory solution if we want 
to make quantum factoring realistic and eventually practical. The authors 
believe that an efficient quantum error correction is needed in order to 
achieve this goal. 
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Figure 1: The Bloch vector s and its projections in the Cartesian coordinate 
system. The field vector BqCz hes in the 2;-direction. When Bq is constant 
the Bloch vector simply rotates around the field vector and its tip follows 
a closed, circular trajectory, which is parallel to the a;?/-plane. If the qubit 
is in an equal superposition of states and 1, then Sz = and the Bloch 
vector itself lies in the a;2^plane. 



Figure 2: Results from a semiclassical stochastic simulation. In (a) and (b) 
the time evolution of Sx and Sy (solid lines) generated by the fiuctuating 
field Bz is shown for two different ensemble members. The dotted lines 
represent the unperturbed evolution. Time t is in units of uJq^. We have 
chosen cq = Ci = 2~^/^ (i.e. Sx = 0). In (c) and (d) we show averages over 
100 and 500 ensemble members, respectively. B^ averages rapidly to zero 
as the ensemble size increases. 



Figure 3: The interaction between qubit and reservoir causes a conditional 
displacement of the field mode. The sign of the displacement depends on the 
qubit state. In the displaced modes the field state oscillates with amplitude 
=t|Ck- This conditional displacement causes a qubit-field entanglement. 
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Figure 4: Decoherence of a single qubit for a one- dimensional field density 
of states for uJc/T = 100 . Time is in units of T~^ and the proportionality 
factor is equal to 0.1. The three decay regimes can be easily identified. 



Figure 5: Decoherence of a single qubit as a function of time and of the 
ratio rj = uJc/T. Time is in units of T~^ and the proportionality factor has 
been set equal to 0.1. Here (a) shows the result of a numerical integration 
of Eq. (|2^ ) for the one-dimensional density of states while (b) shows of the 
exact solution for the three-dimensional density of states (Eq.p4|) 



Figure 6: The decay of two qubit coherence in the case of the shared reser- 
voir with one- dimensional density of states (Eq. ^^. We have set r/ = 1 and 
the proportionality factor is chosen such that at the limit of large ts we get 
the results of Fig. 5(a). In (a) we see how the decay is cancelled out , and 
in (b) we see how it is amplified , when tg is small. The onset of decay does 
not change with ts- 



Figure 7: The decay of two qubit coherence in the case of the shared reser- 
voir with three-dimensional density of states (Eq. ^8|). We have set r^ = 1 
and the proportionality factor is chosen such that at the limit of large t^ we 
get the results of Fig. 5(b). In (a) we see how the decay is cancelled out, 
and in (b) we see how it is amplified, when tg is small. The saturation of 
the decay is present in both cases, though. 
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